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AC  impedance  or  electrochemical  impedance  spectroscopy  (EIS)  is  becoming  a  fundamental  technique 
used  by  researchers  and  scientists  in  proton  exchange  membrane  (PEM)  fuel  cell  analysis  and  develop¬ 
ment.  In  this  work,  in  situ  impedance  measurements  are  presented  for  a  series  of  operating  conditions 
in  a  50  cm2  fuel  cell.  The  electrode  charge  transfer  resistance  was  determined  from  the  corresponding 
arcs  of  the  Nyquist  diagrams.  The  analyses  were  performed  for  H2/02  and  H2/air  operation  at  different 
stoichiometric  factors  and  reactant  gases  humidification.  Characteristic  time  scales  of  charge  transfer 
processes  at  the  different  operating  conditions  were  estimated  from  the  corresponding  Bode  plots.  These 
values  were  used  for  a  non-dimensional  analysis  of  the  different  fuel  cell  electrochemical  and  transport 
processes,  namely  electrochemical  reaction  versus  GDL  reactant  transport.  Fuel  cell  adapted  Damkholer 
numbers  are  thus  presented,  where  the  results  indicate  that  the  GDL  diffusion  transport  is  the  limiting 
process  for  the  cases  under  analysis,  especially  when  air  is  used  as  oxidant.  Additional  analysis  of  channel 
convective  mass  transport  versus  GDL  diffusive  mass  transport  is  also  presented. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

AC  impedance  (electrochemical  impedance  spectroscopy)  is  an 
experimental  technique  widely  used  in  the  analysis  of  electrochem¬ 
ical  systems,  and  has  been  established  as  a  fundamental  diagnostic 
and  research  tool  for  PEM  fuel  cells  during  the  last  years  [1,2].  One  of 
the  main  advantages  of  AC  impedance  techniques  is  that  it  allows 
for  the  determination  of  the  different  resistances  occurring  in  a 
PEM  fuel  cell,  corresponding  to  the  different  electrochemical  and 
transport  processes:  activation  or  charge  transfer  resistance,  ohmic 
resistance,  and  transport  or  concentration  resistance.  On  contrast, 
the  polarization  curve  itself  does  not  provide  detailed  information 
about  the  different  processes  as  they  are  tightly  coupled  and  the 
polarization  curve  describes  only  the  integral  output  of  all  pro¬ 
cesses.  AC  impedance  is  widely  applied  to  multiple  MEA  analysis 
such  as  influences  of  the  catalyst  loading,  PTFE  content,  Nation  con¬ 
tent,  GDL  structure,  manufacturing  methods,  membrane  thickness, 
and  others  [1  ].  In  an  important  number  of  studies  EIS  analysis  has 
been  applied  to  characterize  fuel  cell  processes.  Parthasarathy  et  al. 
[3],  Antoine  et  al.  [4]  and  Neyerlin  et  al.  [5]  analysed  the  oxygen 
reduction  reaction  (ORR)  kinetics.  Xu  et  al.  [6]  and  Neyerlin  et  al. 
[7]  also  analysed  the  effect  of  relative  humidity  conditions  on  the 
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ORR  kinetics  of  PEM  fuel  cells.  HOR  (hydrogen  oxidation  reaction) 
kinetics  and  the  effect  of  CO  poisoning  were  studied  by  Wagner  and 
Giilzow  [8].  Transport  losses  were  analysed  by  Springer  et  al.  [9], 
Eikerling  and  Kornyshev  [10],  Lefebvre  et  al.  [11],  and  Saab  et  al. 
[12,13],  who  also  analysed  the  influence  of  the  processing  con¬ 
ditions  of  the  electrode.  Liu  et  al.  [14]  also  studied  the  resistance 
to  the  protonic  conduction  at  the  cathode  ionomer  layer  using  AC 
impedance.  The  cell  ohmic  resistance  is  also  commonly  determined 
by  means  of  AC  impedance,  as  in  the  work  of  Cooper  and  Smith  [15]. 
Romero-Castanon  et  al.  [16],  Song  et  al.  [17]  or  Wagner  [18]  have 
applied  AC  impedance  for  the  evaluation  and  optimization  of  the 
membrane  electrode  assembly  (MEA).  Additional  applications  are 
the  analysis  of  electrode  and  GDL  flooding  or  membrane  dry-out, 
what  is  commonly  referred  to  as  MEA  State-of-Health  (SOH),  as  in 
the  work  reported  by  Fouquet  et  al.  [19]  or  Merida  et  al.  [20].  The 
results  are  presented  for  the  stack  under  analysis,  but  the  analysis 
of  single  cells  in  a  stack  is  also  feasible  as  reported  by  Hakenjos 
et  al.  [21  ],  who  presented  a  measurement  set-up  featuring  a  multi¬ 
channel  frequency  response  analyser  (FRA)  for  the  simultaneous 
measurement  of  impedance  spectra  of  single  cells  in  a  fuel  cell 
stack. 

Non-dimensional  analysis  has  proven  to  be  a  powerful  tool  for 
engineering  design  and  analysis,  especially  in  the  chemical  engi¬ 
neering  field,  although  its  application  to  complex  systems  such  as 
fuel  cells  is  not  yet  fully  established  [22].  However,  there  is  a  need 
for  comprehensive  and  numerically  less  expensive  description  of 
cell  performance  in  order  to  implement  it  into  full  stack  models 
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Nomenclature 

a  specific  catalyst  surface  area  (m2  m-3 ) 

d  channel  depth  (m) 

D0  diffusion  coefficient  (m2  s-1 ) 

D^pemfc  Damkholer  number 
F  Faraday  constant  (96,485  As  mol-1 ) 

/peak  frequency  (Hz) 

h  mass  transfer  coefficient  (m  s-1 ) 

HFR  high  frequency  resistance  (£2) 

I  current  intensity  (A) 

J  diffusive  flux  (kg  m-2  s-1 ) 

L  channel  total  length  (m) 

Mw  molecular  weight  (g  mol-1) 

P  pressure  (Pa) 

Qv  volume  flow  (mlw  min-1 ) 

R  bulk  resistance  (£2) 

Rc  contact  resistance  (£2) 

Peat  cathode  exchange  current  density  (A  m-3 ) 

So2  oxygen  consumption  (kg  m-3  s-1 ) 

Sh  Sherwood  number 

T  temperature  (I<) 

V  flow  velocity  (ms-1) 

y  mass  fraction 

Greek  symbols 

5gdl  GDL  thickness  (m) 

s  porosity 

p  density  (kg  m-3) 

r  characteristic  time  scale  (s) 

X  stoichiometric  factor 

Subscripts 

0  reference  state 

BP  bipolar  plate 

ch  channel 

chMT  channel  convective  mass  transfer 

chrt  flow  residence  time 

diff  diffusion 

EC  electrochemical 

eff  effective 

GDL  gas  diffusion  layer 

MEA  membrane  electrode  assembly 

MT  mass  transfer 

S  catalyst  surface 


that  cannot  be  analysed  by  means  of  Computational  Fluid  Dynamics 
(CFD)  models  except  for  small-scale  units  [23  ].  Therefore,  a  detailed 
description  based  on  dimensional  analysis  of  the  cell  phenomena 
would  be  of  high  relevance  for  further  stack  research  and  devel¬ 
opment,  and  important  contributions  are  already  in  progress  as 
established  in  the  work  of  Gyenge  [22]. 

This  work  presents  AC  impedance  results  for  a  50  cm2  PEM  fuel 
cell.  The  electrode  activity  and  membrane  protonic  conductivity 
were  determined  elsewhere  [24]  based  on  the  measurements  per¬ 
formed  for  different  operating  conditions,  aimed  at  the  fuel  cell 
parameter  estimation  needed  for  numerical  CFD  simulations  [25]. 
In  this  work,  the  application  of  AC  impedance  measurements  to  the 
non-dimensional  analysis  of  PEM  fuel  cell  phenomena  is  presented, 
and  results  for  a  50  cm2  PEM  fuel  cell  are  reported.  In  particular, 
Nyquist  and  Bode  plots  are  presented  for  a  fuel  cell  with  serpentine 
flow  field,  operating  with  oxygen  and  air  at  different  humidifi¬ 
cation  conditions.  Measurements  at  several  current  densities  are 
presented  in  order  to  better  assess  the  influence  of  the  current 


Fig.  1.  Serpentine  bipolar  plate  used  in  the  analysis. 


density  on  the  limiting  processes  affecting  the  performance  of  the 
cell. 

2.  Cell  description,  operating  conditions  and  experimental 
analysis 

2.1.  PEM  single  cell  description 

The  single  fuel  cell  analysed  consists  of  commercial  graphite 
bipolar  plates  with  a  five-channel  serpentine  flow  field  design 
(Fig.  1)  from  ElectroChem  Inc.  (USA).  The  GDL  is  a  Sigracet  10 
CC  from  SGL  Group  (Germany),  with  420  [xm  thickness,  poros¬ 
ity  0.82  and  10%  PTFE  content,  featuring  a  micro  porous  layer 
(MPL).  The  50  cm2  membrane  used  is  a  CCM  type  (catalyst  coated 
membrane)  from  Baltic  fuel  cells  (Germany),  with  Nafion-117  and 
0.3  mg  Pt/cm2  and  0.6  mg  Pt/cm2  catalyst  load  in  anode  and  cath¬ 
ode  respectively.  The  ratio  of  ionomer/carbon  (I/C  ratio)  is  1.2/2, 
with  70%  Pt  over  carbon,  and  a  ratio  ionomer/catalyst=  1/4. 

2.2.  Operating  conditions 

The  operating  conditions  used  in  the  analysis  are  listed  in 
Table  1.  Cell  temperature  and  pressure  were  maintained  constant 
during  all  the  analysis  at  60  °C  and  4  bar  (absolute),  which  are  the 
operating  conditions  recommended  by  the  cell  supplier. 

2.3.  Experimental  analysis 

The  experimental  work  was  conducted  in  a  FuelCon  CT-1000 
(Germany)  test  station,  located  at  INTA  facilities.  The  experimen¬ 
tal  methods  and  procedures  were  carried  out  following  the  details 
described  in  the  FCTESTNET  documents  for  single  cells  [26],  where 


Table  1 

Operating  conditions  used  for  the  experimental  analysis. 


Aa 

RHa (%) 

Ac 

RHC  (%) 

Oxidant 

1.5 

60 

5.0 

60 

Air 

1.5 

100 

5.0 

100 

Air 

2.0 

100 

10.0 

100 

Oxygen 
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Fig.  2.  Nyquist  diagrams  for  the  operation  at  160  mAcm-2. 


its  particular  application  to  the  case  analysed  is  described  in  [24]. 
AC  impedance  measurements  were  conducted  in  a  frequency  range 
20  Hz  to  20  kHz  on  the  FuelCon  test  station  using  its  own  frequency 
response  analyzer  (FRA).  The  AC  impedance  measurements  were 
conducted  in  galvanostatic  mode,  where  the  measurements  are 
performed  by  applying  an  AC  current  perturbation  to  the  cell  and 
measuring  the  potential  response  of  the  cell. 

It  must  be  mentioned  that  no  measurements  were  carried  out 
in  the  low  frequency  range  (0.1-10  Hz),  although  such  experimen¬ 
tal  work  would  be  needed  in  order  to  obtain  further  information 
about  the  mass  transport  phenomena,  and  also  to  provide  the  nec¬ 
essary  results  for  the  validation  of  the  non-dimensional  analysis 
presented  in  Section  4.  This  will  be  accomplished  in  the  next  stage 
of  the  investigation. 

3.  Results 

The  results  of  the  measurements  are  presented  in  the  Nyquist 
diagrams  shown  in  Figs.  2  and  3.  The  diagrams  feature  two 
impedance  arcs.  The  ZRe  axis  intercept  at  high  frequencies  rep¬ 
resents  the  overall  resistance  of  the  cell  (electronic  at  the  bulk 
materials,  electronic  contact  resistances,  and  ionic  resistances  at 
the  membrane)  and  it  corresponds  to  the  high  frequency  resis¬ 
tance  (HFR).  In  general,  the  next  semicircle  diameter  at  medium 
frequencies  represents  the  activation  overpotential  losses,  and  the 
diameter  of  the  third  semicircle  (at  low  frequencies)  represents  the 
concentration  losses.  However,  the  low  frequency  semicircle  is  typ¬ 
ically  missing  in  cells  operating  with  pure  oxygen  [9]  as  is  the  case 
represented  in  Figs.  2  and  3.  Moreover,  the  measurements  were 
not  conducted  until  the  low  frequencies  at  which  this  semicircle 
would  occur,  as  the  lowest  frequency  reported  in  Figs.  2  and  3  was 
20  Hz,  while  concentration  losses  shall  be  measured  in  the  range 
0.1 -10  Hz.  Therefore,  the  two  semicircles  observed  in  the  Nyquist 
plots  in  Figs.  2  and  3  are  associated  to  the  charge  transfer  resistance 
or  activation  polarization.  The  arc  in  the  left  of  the  diagram  corre- 


Fig.  3.  Nyquist  diagrams  for  the  operation  with  H2/O2  and  full  humidification. 


sponds  to  the  anode  charge  transfer  process  and  the  arc  in  the  right 
of  the  diagram  corresponds  to  the  cathode  charge  transfer  process. 
Although  in  most  cases  both  arcs  overlap  themselves  featuring  a 
unique  arc,  the  results  obtained  show  two  separate  arcs  of  appre¬ 
ciably  different  sizes.  This  indicate  that  the  anode  charge  transfer 
resistance  is  much  smaller  than  the  cathode  resistance  but  still  sig¬ 
nificant,  probably  due  to  the  lower  anode  catalyst  loading.  None 
of  the  arcs  is  representing  the  concentration  polarization  losses  or 
mass  transport  resistance  (the  reduction  of  the  surface  concentra¬ 
tion  of  oxygen,  mainly  caused  by  limitations  in  the  gas  diffusion 
and  pore  blockage  by  liquid-phase  accumulation  at  high  current 
densities  [27]).  This  specific  resistance  would  correspond  to  AC  fre¬ 
quencies  below  10  Hz,  and  measurements  were  performed  in  the 
range  20  Hz  to  20  kHz. 

Fig.  3  shows  the  dependence  of  the  charge  transfer  resistance 
at  different  current  densities  for  the  operation  with  oxygen.  It  is 
clearly  observed  that  the  diameter  of  the  semicircle  associated 
to  the  cathodic  charge  transfer  resistance  dramatically  decreases 
as  the  current  density  increases,  thus  indicating  that  at  low  cur¬ 
rent  densities  the  cathode  activation  losses  are  dominating  the  cell 
voltage  loss.  The  reaction  rate  determining  the  cell  performance 
is  the  kinetics  of  the  oxygen  reduction  reaction  (ORR)  at  the  cath¬ 
ode,  and  thus  the  charge  transfer  resistance  is  determined  by  the 
radius  of  the  corresponding  impedance  arc.  The  relation  between 
the  charge  transfer  resistance  or  activation  polarization  and  the 
electrode  overpotential  is  given  by  the  Tafel  equation.  As  current 
density  increases  the  cathode  overpotential  decreases,  as  there  is 
more  energy  available  for  the  electrochemical  reaction  to  occur.  For 
all  cases  a  significantly  larger  resistance  is  observed  for  the  cath¬ 
ode  rather  than  for  the  anode,  especially  at  low  current  densities,  as 
corresponds  with  the  slower  kinetics  of  the  ORR.  The  reason  for  the 
exponential  decay  in  the  activation  region  of  the  polarization  curve 
is  clearly  observed  in  Fig.  3,  where  the  influence  of  the  current  den¬ 
sity  over  the  cathode  activity  is  shown.  Setting  the  current  density 
from  120 mAcm-2  to  200 mAcm-2  results  in  a  much  larger  elec¬ 
trode  resistance  reduction  (50%)  than  changing  the  current  density 
from  400 mAcm-2  to  500 mAcm-2  (6%). 

As  discussed,  the  electrode  activity  is  determined  from  the  cor¬ 
responding  arc.  Fig.  2  shows  the  dependence  of  the  cathode  charge 
transfer  resistance  on  the  reactant  concentration  and  humidifi¬ 
cation,  where  a  minor  charge  transfer  resistance  is  observed  for 
the  operation  with  pure  oxygen.  When  air  is  used  as  reactant,  an 
enhanced  humidification  increases  the  electrode  activity,  which 
was  also  reported  by  Xu  et  al.  [6]  and  Neyerlin  et  al.  [7]. 

The  total  ohmic  resistivity  of  the  cell  (bulk  electric  resistance, 
electric  contact  resistance  and  protonic  resistance  in  the  mem¬ 
brane)  is  determined  in  the  Nyquist  diagram  as  the  ZRe  value 
at  Zim  =  0  [27,28].  This  value  is  commonly  known  as  high  fre¬ 
quency  resistance  (HFR)  as  it  is  obtained  at  high  AC  frequencies. 
However,  a  single  measurement  at  1  kHz,  which  is  sometimes 
found  as  a  recommendation  in  the  literature  for  the  determina¬ 
tion  of  HFR,  is  not  correct  as  observed  in  Figs.  2  and  3-  Zim  -  0 
is  obtained  at  different  frequencies  always  higher  than  1  kHz  for 
the  cell  under  analysis,  and  the  spectra  show  that  at  1  kHz  the 
impedance  is  not  purely  resistive,  presenting  a  non-zero  imaginary 
component. 

As  the  bulk  electric  resistance  is  known  from  the  materials 
data  sheets,  and  the  BP-GDL  contact  resistance  can  be  deter¬ 
mined  experimentally  [24,28],  the  MEA  resistivity  can  be  calculated 
as: 


^mem  —  HFR  —  RR  P  —  ^GDL  —  2RC  (1 ) 


Results  for  the  cell  under  analysis  were  reported  in  [24]. 
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Frequency  [Hz] 


Fig.  4.  Bode  diagram  corresponding  to  an  operation  with  H2/O2  and  RH  =  100%.  AC 
frequency  range  20  Hz  to  20  kHz. 


4.  Discussion  and  application  to  the  non-dimensional 
analysis  of  fuel  cell  phenomena 

4.1.  Electrochemical  and  GDL  diffusion  mass  transport  time 
scales:  Damkholer  number 


The  results  plotted  in  the  Nyquist  diagram  can  also  be  inter¬ 
preted  as  characteristic  time  scales  of  the  different  processes 
occurring  within  the  fuel  cell.  The  time  constant  r  of  a  particular 
physical  or  chemical  phenomenon  corresponds  to  the  frequency 
peak  of  the  corresponding  arc  in  the  Nyquist  plot.  Therefore, 
r  =  1  //peak  [29].  As  frequency  is  not  explicitly  shown  in  Nyquist  plots, 
the  representation  of  the  impedance  versus  frequency  in  a  Bode 
diagram  is  more  convenient. 

By  determining  the  time  constants  of  the  different  processes  a 
non-dimensional  analysis  can  be  derived.  In  particular,  a  fuel  cell 
Damkholer  number  can  be  used  so  that: 


EC  reaction  rate 
D3pEMFC  = - MTrate - 


(2) 


DapEMFC  =  ~~ 
^EC 


(3) 


The  electrochemical  reaction  rate  (EC  reaction  rate)  corresponds 
to  the  rate  of  the  oxygen  oxidation  reaction  (ORR)  as  it  is  the  limiting 
reaction  in  most  cases  [5,30].  The  mass  transport  rate  (MT  rate)  is 
assumed  to  be  controlled  by  oxygen  reactant  diffusion  in  the  GDL. 
The  derivation  of  the  time  constants  and  the  related  Damkholer 
numbers  are  presented  in  the  following  sections. 


4.1.1.  Electrochemical  time  scales 

Fig.  4  shows  the  Bode  diagrams  for  120  and  500  mA  cm-2  for 
the  cell  operating  with  H2/O2  and  full  humidification,  where  the 
corresponding  Nyquist  plot  was  shown  in  Fig.  3.  The  imaginary 
component  of  the  impedance  is  represented  versus  the  frequency. 
The  peaks  of  the  impedance  values  correspond  to  the  maximum 
values  of  imaginary  impedance  in  the  Nyquist  plots,  and  from  the 
frequencies  featuring  such  peak  impedances  the  time  constants  can 
be  derived.  The  time  constants  for  the  ORR  (dashed  lines)  are  cal¬ 
culated  according  to  r  =  l//peak  [29],  so  that  rEC  can  be  obtained 
from  the  Bode  diagrams.  As  shown  in  Fig.  4,  the  enhanced  driving 
force  for  the  electrochemical  reaction  at  larger  current  densities 
is  involving  a  minor  time  constant:  9.0  ms  for  500  mA  cm-2  and 
25.0  ms  para  120  mA  cm-2. 


4. 1.2.  GDL  diffusion  mass  transport  time  scales 

The  value  of  rdiff  can  be  obtained  in  the  same  manner  from 
Nyquist/Bode  diagrams  representing  the  low  frequency  arc  in  the 
frequency  range  0.1 -10  Hz  [1,19].  Alternatively,  different  approxi¬ 
mations  can  be  used,  such  as  [31  ]: 


^diff  = 


S2 

°GDL 

Deff 


(4) 


where  5gdl  is  the  GDL  thickness  and  Deff  is  the  oxygen  effective 
diffusion  coefficient.  rdiff  represents  the  diffusion  time  or  char¬ 
acteristic  time  scale  of  diffusion,  which  is  the  time  for  significant 
diffusion  to  reach  a  given  distance  8  [32]. 

In  order  to  determine  Deff,  the  standard  diffusion  coefficient 
must  be  corrected  with  the  GDL  porosity  and  the  degree  of  flood¬ 
ing  by  liquid  water  pore  blockage,  where  the  correlation  from  Nam 
and  Kaviany  [33]  has  been  used  in  the  present  analysis.  Addition¬ 
ally,  temperature  and  pressure  effects  must  be  also  considered  as 
stated  in  the  work  by  Um  et  al.  [34].  Considering  all  the  effects 
affecting  the  effective  diffusivity  of  reactants  through  the  GDL,  Deff 
is  calculated  as: 


Vs,^)  (i)' 


(5) 


where  T0,  Po  are  the  reference  temperature  and  pressure,  300  K  and 
101,325  Pa,  and  s  is  the  water  saturation  (volume  fraction  of  liquid 

water).  The  terms  correcting  D0  are  ( fZo’11  )°'785>  correspond¬ 
ing  to  the  correction  for  the  GDL  porosity  and  degree  of  flooding 
(according  to  the  work  of  Nam  and  Kaviany  [33]  for  details),  and 
'  ,  corresponding  to  the  correction  for  temperature  and 
pressure  according  to  Um  et  al.  [34]. 

Another  well  established  model  for  describing  the  effective  dif¬ 
fusivity  in  GDL  porous  media  proposes  that: 


Deff  =  DoS1  5 


(6) 


where  the  exponent  1.5  is  used  in  order  to  account  for  the  medium 
tortuosity.  This  is  known  as  the  Bruggeman  correction  (the  origi¬ 
nal  work  was  published  by  Bruggeman  [35]),  and  the  approach  is 
commonly  used  in  fuel  cell  modelling  [36-38]. 

Both  Eqs.  (5)  and  (6)  yield  similar  results  for  the  high  porosities 
typically  used  in  the  GDL.  For  the  cell  under  analysis  described  in 
Section  2.1 ,  with  s  =  0.82  [39],  differences  in  the  Deff  values  are  less 
than  8%  when  calculated  according  to  Eqs.  (5)  and  (6). 

With  a  GDL  porosity  £  =  0.82,  an  oxygen  diffusivity  in 
air  D0 o2-air  —  1 .88E  —  5  m2  s-1  (average  between  the  value 
reported  in  transport  properties  tables  [40]  and  the  value  result¬ 
ing  from  the  Fuller-Schettler-Giddings  equation),  and  considering 
the  cell  operating  conditions  (4  bar,  60°  C),  the  value  calculated 
using  Eq.  (5)  for  Deff_Air  is  4.6E-6m2s-1,  which  is  25%  of  the 
standard  oxygen  diffusivity  value  in  air.  For  the  case  operat¬ 
ing  with  oxygen,  D0io2-h2o  =  2.27E  -5  m2  s-1  according  to  the 
Fuller-Schettler-Giddings  equation,  and  considering  the  cell  oper¬ 
ating  conditions  the  value  calculated  for  Deff  o2  is  5.56E-6m2  s_1. 
The  water  saturation  value  has  been  considered  to  be  zero  at  the 
current  densities  and  cell  conditions  under  analysis.  This  is  because 
even  at  the  highest  current  density  in  Fig.  3(500  mA  cm-2 ),  no  con¬ 
centration  polarization  effects  are  still  observed  in  the  polarization 
curve  (Fig.  5).  The  calculated  values  of  rdiff  (Eq.  (4))  are  38  ms  for 
air  and  32  ms  for  oxygen. 

For  hydrogen,  D0,h2-h2o  =  8.16E  -  5  m2  s-1  (average 
between  the  value  reported  in  transport  properties  tables  [18]  and 
the  value  resulting  from  the  Fuller-Schettler-Giddings  equation), 
and  considering  the  cell  operating  conditions  (4  bar,  60°  C),  the 
value  calculated  for  Deff_H2  =  2.0E  -  5  m2  s-1 . 
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Fig.  6.  Electrochemical  characteristic  times  for  operation  with  air  and  oxygen. 


Fig.  5.  Polarization  curves  of  the  cell  analysed. 


4.  \ .3.  Damkholer  numbers 

The  characteristic  time  scales  for  the  electrochemical  reaction 
rate  and  the  reactant  diffusion  in  the  GDL,  as  calculated  in  the  pre¬ 
vious  sections,  are  reported  in  Table  2.  The  resulting  Damkholer 
numbers  are  also  presented  in  Fig.  6.  The  results  indicate  that 
the  reactant  diffusion  in  the  GDL  is  the  limiting  process  for  the 
cases  considered,  especially  when  oxygen  is  used  as  reactant  at 
higher  current  densities.  This  result  is  obtained  even  at  low  current 
densities  (160  mA cm-2)  indicating  a  non-balanced  design  of  the 
electrode  and  GDL  for  this  cell,  which  is  probably  a  consequence  of 
the  selection  of  cell  components  from  different  suppliers  (described 
in  Section  2.1). 

The  balance  between  electrochemical  reaction  and  diffusive 
mass  transport  can  be  further  analysed  as  follows.  In  the  steady- 
state  operation,  the  reactant  diffusive  mass  transport  equals  the 
reactant  consumption  via  the  electrochemical  reaction.  The  oxygen 
consumption  at  the  electrode  is: 


balanced  while  one  of  them  is  limiting  the  overall  rate  of  the  pro¬ 
cess.  The  balance  can  be  written  by  combining  Eqs.  (7)  and  (8): 


P^eff-02 

^GDL 


(yo2,GDL 


Vo2,s) 


1  (  Mw,o2 
a  [  4  F 


^cat 


(9) 


For  the  cases  analysed,  the  limiting  term  in  Eq.  (9)  is  the  GDL 
diffusive  mass  transport  as  demonstrated  by  the  corresponding 
Damkholer  numbers. 


4.2.  Channel  and  GDL  mass  transport  time  scales 

Channel  design  must  ensure  an  appropriate  reactant  distribu¬ 
tion  over  the  electrode.  Given  the  typical  channel  dimensions  and 
reactants  mass  flows,  the  flow  regime  is  always  laminar,  so  that  the 
convective  mass  transport  is  represented  by  the  Sherwood  number 
[41]: 

Shch  =  ^~2.6  (10) 


So2=- 


f  ATw,o2 

^  4  F 


Feat 


(7) 


where  Rc at  is  the  cathode  volumetric  transfer  current  given  by  the 
Butler-Volmer  equation,  Mw?c>2  is  the  oxygen  molecular  weight, 
and  F  is  the  Faraday  constant.  The  diffusive  flux  from  the  GDL  to 
the  electrochemical  active  surface  area  is  defined  by: 

Jo2  =  —  /ff  °2  (y02,GDL  —  yo2,s)  (8) 

aGDL 

where  p  is  the  density,  y02,GDL  an^  Yo2,s  is  th e  oxygen  mass  fraction 
at  GDL  and  reaction  sites  respectively.  The  volumetric  consumption 
of  oxygen  So2  can  be  transformed  into  a  surface  flux  by  dividing 
its  value  by  the  specific  surface  area  of  the  catalyst  or  surface- 
to-volume  ratio  (a).  In  steady  state  operation,  both  processes  are 


where  h  is  the  convective  mass  transfer  coefficient  for  channel- 
GDL,  d  is  the  channel  depth,  and  Dch  is  the  diffusion  coefficient.  The 
characteristic  time  scale  for  the  reactant  convective  mass  transport 
is  therefore  d/h,  or: 

d  d2 

^chMT-/j~2 ,6Dch  *■  ^ 

For  the  channel  depth  of  the  bipolar  plate  shown  in  Fig.  1 
(d  =  0.001 1  m),  rchMT  =  20  ms  for  02, 25  ms  for  air,  and  57  ms  for  H2. 

By  combining  Eqs.  (4)  and  (1 1 )  a  characteristic  time  scales  ratio 
for  the  reactant  transport  in  channels  (convective)  and  GDL  (diffu¬ 
sive)  can  be  obtained: 

rchMT  =  ^GDL 

tdiff  2.6Dc[1Deff 


Table  2 

Electrochemical  reaction  characteristic  time  scale,  Damkholer  numbers,  and  flow  residence  times.  5GDL  =  0.42mm,  Deff.air  =  4.6E-6m2  s_1,  Deff_o2  =5.56E-6  m2  s-1, 
Tdiff-air  =  38  ms,  Ddiff_o2  =  32  ms,  rchMT  =  20  ms  for  02,  25  ms  for  air,  and  57  ms  for  H2. 


Aa 

RHa  (%) 

Ac 

RHC  (%) 

Oxidant 

i  (mAcur2) 

rEc  (s) 

DaPEMFC 

"tchrt.  3  (s) 

tchrt.  C  (s) 

1.5 

60 

5.0 

60 

Air 

160 

0.033 

0.87 

1.939 

0.243 

1.5 

100 

5.0 

100 

Air 

160 

0.028 

0.74 

1.939 

0.243 

2.0 

100 

10.0 

100 

Oxygen 

120 

0.025 

0.75 

1.939 

0.774 

2.0 

100 

10.0 

100 

Oxygen 

160 

0.020 

0.63 

1.454 

0.581 

2.0 

100 

10.0 

100 

Oxygen 

200 

0.016 

0.52 

1.163 

0.465 

2.0 

100 

10.0 

100 

Oxygen 

400 

0.010 

0.21 

0.582 

0.232 

2.0 

100 

10.0 

100 

Oxygen 

500 

0.009 

0.28 

0.465 

0.186 
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A  combined  design  of  both  channel  and  GDL  should  ensure 
that  rchMT/rdiff  ~  1,  so  that  no  limiting  process  is  controlling  the 
mass  transfer  rate.  However,  it  is  widely  known  that  in  general 
TchMT  <  rdiff  given  the  GDL  characteristics  in  terms  of  porosity  and 
water  management,  especially  at  high  current  densities.  For  the 
cases  considered,  rchMT/rdiff  =  0.65,  indicating  a  mass  transfer  pro¬ 
cess  controlled  by  the  diffusive  mass  transport  in  the  GDL. 

In  the  same  manner,  the  flow  residence  time  in  the  channel  must 
ensure  that  the  convection  mass  transport  occurs,  thus  the  flow 
residence  time  must  be  larger  than  the  convective  mass  transport 
time  scale: 

^chrt  >  ^chMT  (13) 

where  rChrt  =  L/V.  For  the  bipolar  plate  considered,  L  =  0.63  m. 
The  flow  velocity  V  depends  on  the  channel  cross-section 
(0.78  mm  x  1.1  mm)  and  the  volume  flow,  which  depends  on  the 
current  density  and  stoichiometric  factor: 

Qv,h2  =  6.97/Aa  (14) 

Qv,o2  =  3.49/Ac  (15) 

0v,Air  —  16.7/AC  (16) 

where  I  is  the  current  intensity  (A)  and  X  the  stoichiometric  fac¬ 
tor.  For  the  cell  under  analysis  with  five  path-serpentine  flow  field 
design,  the  total  Qy  is  split  into  the  five  parallel  channels  in  order 
to  analyse  a  single  channel.  Table  2  lists  the  rchrt  values  obtained 
for  this  analysis  for  the  anode  and  cathode  sides.  The  comparison 
between  rchrt  values  against  rchMT  values  shows  that  the  residence 
times  of  anode  and  cathode  flows  are  long  enough  in  order  to  assure 
the  appropriate  mass  transfer  process  to  the  GDL  for  all  cases  con¬ 
sidered  in  the  analysis. 

5.  Conclusions 

AC  impedance  results  for  a  50  cm2  PEM  fuel  cell  are  presented. 
Besides  analyzing  the  electrode  charge  transfer  resistance  from  the 
Nyquist  plots  obtained  for  different  operating  conditions,  this  work 
outlines  a  non-dimensional  analysis  of  PEM  fuel  cell  transport  phe¬ 
nomena.  Characteristic  time  scales  for  electrochemical  reaction  and 
reactant  transport  processes  in  channels  and  GDL  are  derived  either 
from  experimental  Bode  plots  or  from  correlations  found  in  the 
literature.  Damkholer  numbers  were  calculated  for  a  range  of  oper¬ 
ating  conditions,  which  can  help  understanding  and  quantifying 
the  fundamental  processes  limiting  the  performance  of  a  PEM  fuel 
cell.  Characteristic  times  for  channel-GDL  convective  mass  transfer 
and  GDL  bulk  diffusive  mass  transfer  was  also  compared,  and  the 
appropriateness  of  the  channel  design  in  terms  of  flow  residence 
time  was  studied.  For  the  cell  and  operating  conditions  analysed, 
the  investigation  shows  that  diffusion  in  GDL  is  constantly  the  lim¬ 
iting  transport  process  even  at  low  current  densities  ( 1 60  mA  cm-2 ) 
indicating  a  non-balanced  design  of  the  electrode  and  GDL  for  this 
cell,  which  is  probably  a  consequence  of  the  selection  of  cell  com¬ 
ponents  from  different  suppliers.  Additional  measurements  in  the 
low  frequency  range  (0.1 -10  Hz)  are  needed  in  order  to  allow  for 
additional  experimental  analysis,  which  will  provide  further  infor¬ 
mation  about  the  mass  transport  phenomena  also  necessary  for  the 
validation  of  the  non-dimensional  analysis  presented  in  this  work. 
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